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Abstract 

We perform a global fit of data relevant to eeqq contact interactions, including 
deep inelastic scattering at high from ZEUS and HI, atomic physics parity 
violation in Cesium from JILA, polarized e~ on nuclei scattering experiments 
at SLAC, Mainz and Bates, Drell-Yan production at the Tevatron, the total 
hadronic cross section cihad at LEP, and neutrino-nucleon scattering from 
CCFR. With only the new HERA data, the presence of contact interactions 
improves the fit compared to the Standard Model. When other data sets are 
included, the size of the contact contributions is reduced and the overall fit 
represents no real improvement over the Standard Model. 



I. INTRODUCTION 



The reports of event rates above Standard Model (SM) expectations in e~^p e^X deep 
inelastic scattering (DIS) at very high by the HI and ZEUS p[ experiments at HERA 
have generated a considerable amount of theoretical activity. Various models to explain the 
excess events have been put forward [^0- Particular attention has been given to resonant 
production of leptoquarks , squarks with _R-parity violating couplings , quark and 
lepton compositeness ||^, and a general parametrization in terms of eeqq contact interactions 
PIJ^JgHTSl, representing the exchange of very massive particles in the s, t or -u-channels. 

An interest in leptoquarks derives from an indication for an e^j invariant mass peak 
at Mej ~ 200 GeV in the HI data. However, no clear mass peak is present in the ZEUS 
data; see discussion in Ref. |T^. In addition, for an ej branching fraction of unity, recent 



searches at the Tevatron exclude leptoquarks of mass Mlq < 210 GeV from CDF and 
Mlq < 225 GeV from DO |TD|. The Tevatron bounds can be escaped, however, with a 
reduced branching fraction B{LQ —>■ ej). 

Any new physics in eg —>■ eq scattering for which the exchanged particles have mass 
squared ^ s can be described by an effective eeqq contact term Lagrangian. For 
example, effects of a Z' boson of TeV mass scale would be well represented by a four- 
fermion contact interaction. A leptoquark of mass M^q ^ 200 GeV could similarly be 
described at HERA by an effective eeqq contact term. For leptoquark masses not far above 
the HERA energy range additional form-factor effects would have to be considered, but an 
analysis in terms of contact interactions would still be appropriate to derive constraints from 
low energy experiments. 

In this paper we perform a combined analysis of low and high energy data to obtain 
constraints on eeqq contact interactions. The data sets include deep inelastic scattering at 
high from ZEUS and HI atomic physics parity violation polarized e~ on nuclei 
scattering experiments at SLAG Mainz |2^, and Bates p3|, Drell-Yan production at 



the Tevatron , the total hadronic cross section (Tj^ad at LEPl, LEP1.5, and LEP2 |2|^, 



the left-right asymmetry Alr at SLD |25], and neutrino-nucleon scattering from GGFR [30 

The organization of the paper is as follows. In Section II we discuss the parametrization 
of contact interactions and discuss possible symmetry relations among them. In Section HI 
we present all the data that enter our global analysis. In Section IV we give constraints on 
the eeqq contact interactions in selected models and show how each experiment constrains 
different combinations of these interactions. Allowing for all eight eeuu and eedd contact 
interactions to vary freely, the high HERA data and the high mass Drell-Yan pair data 
from GDF cannot be fitted simultaneously. Section V gives our conclusions. 



II. PARAMETRIZATION OF CONTACT INTERACTIONS 

The conventional effective Lagrangian of eeqq contact interactions has the form [pT1 -|33 
Lnc = Y.[^LL {elliieL) {qLl^QL) + Vrr [e^'^^^en) {W1^<1r) 

+ Vlr (eZT/^ei) {qRl^qR) + r]RL (e^T^ei?) (QlI'^Ql) , (1) 



2 



where the eight independent coefficients 77^^ and r]^p have dimension (TeV)~^ and are con- 
ventionally expressed as rj'^p = eg'^/A^^, with a fixed = in. The sign factor e = ±1 allows 
for either constructive or destructive interference with the SM 7 and Z exchange amplitudes 
and represents the mass scale of the exchanged new particles, with coupling strength 
g^/4:7[ = 1. A coupling of this order is expected in substructure models and therefore Aeq 
is sometimes called the "compositeness scale" 0]. In the effective interaction (|l|) we do 
not include lepton or quark chirality violating terms such as (eTe/j) {JilQr) since, if there 
is an approximate SU{2) x U{1) invariance, scalar and tensor contact terms are severely 
constrained by helicity suppressed processes like vr" — > e~u A contact term Lagrangian 
completely analogous to Eq. (|1|) can be written down for uugq interactions, with parameters 
?7^^ and r^^^, respectively. 



A. SU{2) X [/(I) symmetry and universality 



Left-handed electrons and quarks belong to SU(2) doublets L = {ul, cl) and Q = {ul, di) 
and thus from SU(2) symmetry one expects relations between contact terms involving left- 
handed u or d quarks; similarly, contact terms for left-handed electrons and neutrinos should 
be related [|14|. In order to exhibit these relationships we start from the most general 
SU{2) X f/(l) invariant contact term Lagrangian, 

+ ??7(e^7^e/j)(c/ij7^c/K) . (2) 

Here we have suppressed the lepton chirality violating term (^L^y'^Q^ {dR'^^^eiij which, under 

a Fierz transformation, is equivalent to a (^Le^j [d^Q^ term which, as noted before, is 
severely constrained by tt^ — e~z/ data. In the second term of Eq. = (t"/2 denotes 

the SU(2) generators; this term describes the exchange of isospin triplet quanta between the 
lepton and quark fields. Because of this possible contribution no constraints arise for the 
contact terms and the difference r/l^, — rjf^i^ = 772/2 measures the size of evud contact 
terms in CC processes. SU(2) does constrain right-handed electron couplings, 

Vrl = V5 = Vrl ■ (3) 
In addition, the four neutrino and the lepton couplings are related by SU(2), 



vll 


= Vi 


1 

+ 


= vfL 


vll 


= Vi 


1 


= Vll 


vTr 


= V3 


= vTr , 




vI'r 


= Va 


= vI'r ■ 





In our analyses, the relations of Eqs. (|^) and (|[) are only used when neutrino scattering 
data are included in the analysis. Even though we expect that SU{2) x f/(l) will be a 
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symmetry of the renormalizable interactions which ultimately manifest themselves as the 
contact terms of Eq. (Q), electroweak symmetry breaking may break the degeneracy of SU(2) 
multiplets of new, heavy quanta whose exchanges give rise to (|1|). This would result in a 
violation of the relations of Eqs. (^) and @. One example is the exchange of the stop ti, 
t2, and the sbottom 6L,&i? in -R-parity violating SUSY models. The large top-quark mass 
may lead to substantial splitting of the masses of these squarks which could easily lead to 
violations by up to a factor of two of SU(2) relations such as 77^"^ = rjf^n- 

In the discussion above we have considered first generation quarks and leptons only, 
because the "HERA anomaly" raises particular interest in such couplings. In principle, all 
r/'s carry four independent generation indices and may give rise to other flavor conserving and 
flavor changing transitions. Because of severe experimental constraints on intragenerational 
transitions like K ^ we restrict our discussion to first generation contact terms. 



Only where required by particular data (e.g. the LEP data which sum over the quarks of all 
three generations) will we assume universality of contact terms between e, and r and up- 
and down-type quarks of different generations. 

Other symmetry constraints on the parametrization of four-fermion contact interactions 
are of interest as well: 



i) SU(12) symmetry considerations |TI| give 



VaL = -VaR , a = L oi R . (5) 

In addition, S'U{2)l gauge invariance requires the w-quark contribution equal to the 
d-quark contribution, i.e., 

vTl = <i , a = L or R. (6) 
With these symmetries the atomic parity violation constraint is satisfied naturally. 



ii) Vector-vector (VV) interactions [12|, from the exchange of a boson with purely vector 
couplings, give 

eq eq eq en eq /rr\ 

Vll = Vrr = Vlr = Vrl = Vvv ■ iV 



iii) Axial-axial interactions |T2|,|T3[, from the exchange of purely axial- vector couplings, 
give 

eq eq eq eq eq /o\ 

Vll = Vrr = -Vlr = 'Vrl = Vaa ■ (8) 

We consider fits with these restrictions imposed subsequently to make model-independent 
analysis. 
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B. Observables and four-fermion amplitudes 



All the observables that we consider are described by a four-fermion ^-matrix element. 
The amplitudes for observables such as e+e~ hadrons, pp —>■ i~^i~X and atomic parity 
violation, are obtained from the amplitude for eq —>■ eq scattering by crossing; they are given 
by angular factors times reduced amplitudes M^^, where the subscripts label the chiralities 
of the initial lepton (a) and quark The SM tree level reduced amplitude for eq — >• eq is 

M^ii) = + ■ , a,(3 = L,R (9) 

^ t sm^ 0^ cos^ 0„t - rriz 

where i = —Q'^ is the Mandelstam variable, = T^f — sm^OvjQf and = —sin^O^Qf, 
Tsf and Qf are, respectively, the third component of the SU(2) isospin and the electric 
charge of the fermion / in units of the proton charge, and = ATcaem- For e+e~ —>■ qq or 
qq — > e~^e~ , the amplitude M^^(s) is obtained from Eq. (H) with the replacement t ^ s and 
t — M| — * s — M| + isTz/Mz, where s is the subprocess cm. energy squared. The new 
physics contributions to the reduced amplitudes Ma/3 from the eeqq contact interactions of 
Eq. (D are 

AM^^ = , a,P = L,R. (10) 



C. Contributions to Standard Model parameters 



In low energy neutral current (NC) processes, at ^/s <^ m|, Z boson exchange can also be 
described by effective four-fermion contact terms. For the parity violating contributions to 
eeqq NC interactions the NC Lagrangian is conventionally expressed in terms of parameters 
Cig and as 



e Hadron 



G 



(11) 



The radiatively corrected SM values are 



^2q 



eq 



-Tsq + 2Qg{K sin^ 9 



~T3q Peq 



1 — 4 sin^ 9^ 



eq 
eq '• 



where sin 9^ 



0.9884 and 



eq 



+ X2q , (12) 

1.036 for the atomic parity violation 



0.2236, p'^g 

experiments while p'^^ = 0.979, k'^^ = 1.034, p^g = 1.002, K^g = 1.06, = —0.013, and 
-^2^ = 0.003 for the lepton-nucleon scattering. The contact term Lagrangian introduces 
parameter shifts 



AC 



AC. 



2V2GF 
1 



eq , eq eq , eq 

Vll + Vrr - Vlr + Vrl 



2q 



2V2G1 



eg I eg I eg 

- Vll + Vrr + Vlr 



eq 

Vrl 



u, d 



(13) 
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that can be constrained by experimental data on atomic parity violation and electron nucleon 
scattering. 

If SU{2) X f/(l) is a good symmetry of the contact interactions then neutrino-nucleon 
scattering data also constrain the eeqq contact terms [Q. Comparing to the conventional 
parametrization of neutrino-quark effective interactions, 



G r 

^ = ^^^C^ - 1^)^ dl 57^(1 - 75)g + 9r 97^(1 + 75)g 

y/ Z q I 

the contact interactions (|l|) introduce shifts in the coefficients g1 ^, 

1 . „ 1 



(14) 



'2V2G, 



eq 

-Vll 



R 



'2V2G, 



eq 

-Vlr 



q = u,d . 



(15) 



Here we have used the SU(2) x f/(l) relations of Eq. (^). The SM values for these couplings 
are 



Wl) 



1.0095 (Tsg - (1.0382 sin^ Q^) + A,l 
{g^f^ = 1.0095 (-g, (1.0382 sin2e^)+Ag^j) , 

with A„L = -0.0032, A^l = -0.0026 and Kr = XdR/2 = 3.6 x 10"^ 



(16) 
(17) 



D. Connection to specific models 

A second Z-boson, Z2, of mass and with chiral couplings eg^^'^^ to fermion species 
/, would give rise to four-fermion contact terms with 

e(2) 9(2) 
eq _ U (3 

Va/3- e , (iSj 

in experiments at s <^ Mf^ . For subprocess energies of the order of the Z2 mass, as may be 
achievable at the Tevatron, form-factor effects, due to the s dependence of the Z2 propagator, 
would have to be considered. 

The exchange of scalar leptoquarks would give rise to contact terms 

eu ^ ^2 ( hlL^aL^m h^RSgRSfSL \ ed ^ ^2 ( hl^^^gL^pR /^ij;(5afl^/3L \ /.gx 

\2{-MlJ^ 2{-Mlj) ' \2{-Miy 2(-Mlj] ' ^''^ 

where eh2L,eh2R,eh2L are the couplings of R2L, R2R, R2L, respectively [^. A leptoquark 



mass in the 200-300 GeV range would again necessitate the inclusion of form-factor (propa- 
gator) effects in the analysis of HERA and Tevatron data. For the low-energy data, however, 
the contact term parametrization would be entirely adequate. 



6 



III. EXPERIMENTAL DATA SETS 



A. ZEUS and HI: e+p e+X 



The HI and ZEUS experiments have provided event rate measurements in x, y bins where 
X and y are the usual scaling variables of DIS: 



X 



QV [2p ■ (A: - A;')] , y = QV[sx]. 



(20) 



Here k and k' are the four-momenta of the incoming electron and outgoing positron, respec- 
tively, p is the proton four-momentum, s ~ (300 GeV)^ is the square of the cm. energy, and 
is the square of the momentum transfer 



{k — k') = sxy . 



(21) 



The observed distribution of events in x, y bins is given in Table I for ZEUS and in Table II 
for HI along with SM expectations that include experimental efficiencies. Systematic 
errors are correlated between bins. However, these correlations have not been published 
yet and we assume their effects to be small compared to the large statistical errors of the 
published data. 

Neutral current deep inelastic scattering occurs via the subprocess eg eg. In terms 
of the reduced amplitudes M^'plJ:) {a, (3 = L,R) of Eq. (||) the spin- and color-averaged 
amplitude-squared for e~^{k)q{p) e^{k')q{p') is given by 



e g 



(22) 



where t = —sxy, u = —sx{l — y) and s = sx. In our analysis we calculate the SM tree 
level cross section for each x, y bin. We then compare with the efficiency corrected SM event 
rates, N^^{x,y; corrected), in Tables I and II to obtain correction factors for each bin. 



C{x,y) 



N {x, y; corrected) 



SM/ 



x,y) 



(23) 



Then in analyzing models that include new physics contributions we multiply the theoret- 
ical cross sections by C {x, y) to take into account the experimental efficiencies and compare 
the corrected result with the observed number of events. We treat the HI and ZEUS data 
separately. The SM DIS differential cross section is given by 



da{e^p) sx 
dx dy 



167r 



MTRii) + M^^^ii) +(1-^)2 M-(t) + (t) 



+ d{x,Q^ 
+m(x,Q2) 

+d{x, Q2) 
+s{x,Q'^) 



Mint) + Ml,l{t) +{l~yyiMll{t) + M-(t) 



MlUi) ' + Mf^(i) \il-yr( Mini) ' + Ml^^ii) 



+c(x,Q2) 



Mrnit) + (t) +il-yy^MlUt) + M^t) 
Mini) ' + Mf^(t) \{l-yf( Mrnii) ' + (t) 



where Q^), di^x^Q"^) etc. are parton distributions. 



(24) 



B. Atomic Parity Violation Experiment 

Parity violation in the SM is due to exchanges of weak gauge bosons, with vector- 
axial-vector {VA) and axial-vector-vector {AV) terms contributing. Atomic physics parity 
violation measures |^ electron-quark couplings that are different than those probed by high 
energy experiments and thereby provides alternative constraints on new physics. In terms 
of the coefficients Ciq and of the NC Lagrangian (|TT]), the weak charge Qw of a heavy 
atom is given by [^| 



Q 



w 



CU2Z + N) + Cid{Z + 2N) 



(25) 



where Z and are the number of protons and neutrons respectively in the nucleus of the 
atom. Recently a very precise measurement was made of the parity violation transition 



between the 6S and 7S states of Cesium with the use of a spin-polarized atomic beam [20 



For ^iCs, 
be 



-376Ci„ — A22Cid- The weak charge Qw of the atom was determined to 



Q^P = -72.11 ±0.27 ±0.89 , (26) 
where the first uncertainty is experimental and the second is theoretical. The new measure- 



ment of (26) represents substantial improvement from the value given in the 1996 Particle 
Data Book [^] and shows better agreement with the SM value. Including radiative correc- 
tions the predicted SM value for Cs is 



qSM ^ -73.11 ±0.05. 
We find the following constraint for the new physics contribution 

AQv^(Cs) = -376ACi„ - 422ACid = 1.00 ± 0.93 . 



(27) 



(28) 



Some chirality combinations of LL, RR, LR, RL give zero contributions to ACig and thus 
satisfy the experimental Qw constraint. Such possibilities include (i) LL = RR = LR = 
RL (VV), (ii) LL = RR = -LR = -RL (AA), (iii) LL = -LR,RL = -RR (an SU(12) 
symmetry ||lll), (iv) LR = RL, LL = RR = (a minimal choice used in fitting the HERA 
data lOl). 
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C. Polarization Asymmetries in Electron-Nucleus Scattering Experiment 



In this subsection, we review three experiments on polarized electron-nucleus scattering: 
the SLAG e-D scattering experiment the Mainz e-Be scattering experiment ||2^, and 



the Bates e-C scattering experiment ||23 



1. SLAC e-D experiment 
The SLAC polarized e-D experiment measured the parity-violating asymmetry 



A 



dan + dai 



(29) 



where da^/L are the differential cross sections for e^^jAD eX scattering. At low energy A 
is given in the valence approximation to the parton model by p5|j37[] 



3Gf 



A 



5V2: 



Ci„- — 



c, 



2u 



C2d 



(30) 



where > is the momentum transfer, y is the fractional energy transfer from the electron 
to hadrons, and Cig, are the coefficients of the parity- violating Lagrangian in Eq. (^TJ). 
After taking into account uncertainties in the sea-quark contribution, the R = ctl/ (Jr ratio 



and the higher-twist effects, the following constraints are found |39 
2Ci„ - Ci, = -0.94 ± 0.26 , 2C2„ - C^d = 0.66 ± 1.23 , p„ 



-0.975, 



(31) 



where pcorr is the correlation. The recent update of the SM contribution to the effective 
parameters at {Q"^) ~ 1.5 GeV^ finds 



SM 



(2Ci„ — Cid) 
(26*211 ~ C2d)^^ 



-0.723 ±0.005, 
0.105 ±0.006. 



We hence use the following constraints for the new physics contributions 
A(2Ci^„-Ci^,) = -0.22 ±0.26, A{2C',^ - C',,) = 0.77 ± 1.23 , p,,„ 
Note the strong negative correlation of the errors. 



(32a) 
(32b) 



-0.975. (33) 



2. Mainz e-Be experiment 

At Mainz, asymmetry of quasi-elastic polarized electron scattering on a ^Be target 
has been measured The experiment was performed at the mean momentum transfer 
(Q^) ~ 0.2025 GeV^. The measured asymmetry parameter ^Mainz is given by the following 
combination of the coefficients of the effective Lagrangian (|n]) 

v4Mainz = 2.73Ci„ - 0.65Cid ± 2.19C2. - 2.03C2d, (34) 
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and the experimental result is 

^Mainz = -0.94 ±0.19. (35) 

Here the error is obtained by taking the quadratic sum of the theoretical and experimental 
errors as quoted in Ref. p2| . 

The theoretical prediction in the SM has been evaluated in 153] to be 



^M'finz = -0.875 ±0.014. (36) 

We hence find the following constraint on new physics contribution 

A^Mainz = 2.73ACi^„ - 0.65ACi^^ ± 2.19ACL - 2.03AC2^, 

= -0.065 ±0.19. (37) 



3. Bates e-C experiment 

The asymmetry of elastic polarized electrons scattering on a ^^C target has been mea- 
sured at Bates . Their asymmetry parameter can be expressed in terms of the effective 
Lagrangian coefficients as follows; 

^Bates=^^[l^^,^(C^l«±C'i,). (38) 

The typical momentum transfer of the experiment is small, (Q^) = 0.0225 GeV^, and con- 
tributions from quarks other than the u and d quarks can be safely neglected ||4^. From the 
measurement 

^Bates = (1.62 ±0.38) X 10^^ (39) 

and the estimate |41] l/a(— 0.0225 GeV^) = 135.87, we obtain the following constraint on 
the effective Lagrangian coefficients, 

<^i«±Cid = 0.137 ±0.033. (40) 

The theoretical prediction of the SM at (Q^) = 0.0225 GeV^ is [0 

{Cm + Cidf^ = 0.1522 ± 0.0004. (41) 
Therefore the constraint on new physics is 

A(Ci'„ ± C^a) = -0.0152 ± 0.033. (42) 
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D. CDF Drell-Yan Cross Sections 



The CDF experiment at the Tevatron has reported prehminary data p3 on Drell-Yan 
lepton-pair production. The data that we use in our analysis are extracted from a CDF 
figure; the values are given in Table |T|. The data are differential cross sections versus M, 
integrated over |?/| < 1 and divided by two to average over rapidity, where M and y are, 
respectively, the invariant mass and rapidity of the lepton pair. The double differential cross 
section versus M and y of the lepton pair is given by 



RR^ 



+ \Ml%{s)\' 



+ \Ml\{s)\' , (43) 



where q{x) and q{x) are parton distributions evaluated at = M^, xi^2 = Me^''^/y/s, 
^/s = 1800 GeV, and s = M^. The reduced amplitudes M^^(s) are given by Eq. (H) with t 
replaced by s. The QCD i^-factor of Drell-Yan production is given by p2[ 



K = 1 + 



a.(M2) 4 



27r 



(44) 



With this K factor, the overall cross section normalization agrees with the CDF data in the 
vicinity of the Z-peak. 



E. e~^e Experiments 



The eeqq contact interactions probed at HERA and in low energy experiments involve 
only light quarks. In e^e~ — >■ hadrons both light and heavy quarks contribute and separation 
of the light quark contributions requires bottom and charm quark tagging. To investigate 
the influence of contact terms on e^e~ hadrons we assume here that the eeqq contact 
interaction is flavor independent. The relevant measurements are the total hadronic cross 
sections Uhad and the left-right asymmetry Alr at LEPl and SLD and ahad at LEP1.5 and 
LEP2. Our analysis is based on the data in Refs. p5|- p9|] . 

At leading order in the electroweak interactions the total hadronic cross section summed 
over all flavors q = u,d, s,c,b is 



(45) 



where the QCD K factor is given by K 
left-right asymmetry A^/j is given by 



l + a,/7r + 1.409(a,/7r)^-12.77(a,/7r)^ ||. The 



A 



LR 



ctl + CTr 



E,(iMr,(.)p + iMMp- 




Mr'r{s)\') 








Mn^RisW) 



(46) 



At ^/s = Mz, the SM amplitudes are imaginary whereas contact term contributions are 
real. Because of the absence of interference with the SM amplitudes, the Z-pole data have 
little sensitivity to the contact terms despite their high accuracy. To take into account 
next-to-leading order (NLO) electroweak radiative corrections we calculate 
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theory _ theory / "NLO | (a^\ 



where (7^0°^^^ is given by (^) with contact interaction contributions included in the matrix 
elements. Table [V| gives the measured (Jhad and the calculated values of cr^i^o- The radiative 
corrections to the LR asymmetry cancel in the ratio (^61) so here we compare the observed 
Alr with the leading order SM calculation. The experimental and the SM predicted values 



for Alji are p5 



A^'^P = 0.1542 ± 0.0037 , A^^^ = 0.145 ± 0.001 ± 0.001 . (48) 

Other LEP1.5 and LEP2 measurements are summarized in Table The new physics 
contributions are given by Eq. ([T0|). 



F. Neutrino-Nucleon DIS Experiments 

Deep inelastic scattering (DIS) experiments with neutrino and anti-neutrino beams have 



provided important tests for the SM since the early 80's p4| , ^ . If SU(2)i invariance is 



assumed, then uN DIS data also constrain contact interactions. The CCFR collaboration 



obtained a model-independent constraint on the effective uugq couplings fsy] 



K = 0.5820 ± 0.0041 = 1.7897^2 + l.U79gl - 0.0916^1 - 0.0782^ , (49) 

where gl j^ = [gl^uf + (^^,/^)^ ^l,r = {qI^rT - {di^Rf and gl,g%^ are the coefficients of 
the effective Lagrangian (14). The Standard Model value is k = ksm = 0.5817 ± 0.0013 
30| . This model- independent constraint ( |i9| ) can be used to constrain physics beyond the 
SM. Since this CCFR result is so far the most accurate measurement for vN scattering, we 
adopt their result in our analysis. We can write g^ = {gDsM + ^di^ ^^'^■) then k becomes 

K = KSM + 1.7897 {2{gl)sMAgl + (Agl)' + 2{gi)sMAgi + (Agi)') 
+ 1.1479 (2((7]^)smA(7]^ + {Ag-^)' + 2(<?^)smA<7^ + {Ag'^)') 

- 0.0916 {2{gl)sMAgl + {Aglf - 2{g'[)sMAg'[ - {Aglf) 

- 0.0782 {2{gl)sMAg-R + {Aglf - 2{g%)^MAgi - {Agif) , (50) 



where the Ag^''^, are the contact term contributions of Eq. (|T^). 



G. Statistical Analysis 



The method of maximum likelihood is the most general method of parameter estimation. 
For a set of independently measured quantities Xj that come from a probability density 
function f[x] a), where a is a set of unknown parameters, the method of maximum likelihood 
consists of finding the set of a, which maximizes the joint probability density for all data, 
given by the likelihood function C = Hi f{xi;a). Very often it is easier to maximize the 
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logarithm of this hkehhood function by solving the equation d\ogC/da = 0. In the case of 
ZEUS and HI data, since the number of observed events in each bin is small, it is appropriate 
to describe the probabilities of the observed events using Poisson statistics: 

_ (nf ) e 

"'i ■ 

where nf^ is the expected number of events from the theory and nf^^ is the number of 
observed events in the ith bin. With Poisson statistics the method of maximum log- likelihood 
is equivalent to minimizing the following 

= E - + K'^ log ■ (52) 

In the bins where nf'^ = the second term is zero. In our fitting procedures, we shall 
vary unknown n parameters of the contact interaction and calculate the expected number 
of events, nf^, for each bin i; then we obtain the quantity fo^' relevant ranges of the 
unknown parameters. The set of parameters that gives the minimum value is the best 
estimate of the parameters. Relative merits of different new physics interactions can also be 
compared by their values. 

In experiments with good statistics, we use the usual approach to calculate 

^ ^ , (53) 

where x1^^, xf^, and a^^ are the experimental measurements, theoretical predictions, and 
the errors, respectively. The method of maximum log likelihood is equivalent to the usual 
method when the distribution follows Gaussian statistics. The x^'s from all data sets 
are added together to form the total x^ and best estimate a for the set of parameters is 
obtained by minimizing the total x^? i-e., X^(^) = Xmm- The n-sigma standard deviation of 



the best estimate a is calculated by finding a', which satisfies 

x\&')=xLn + n^ ■ (54) 

For a l(2)-parameter fit a 1-a standard deviation corresponds to a confidence level of 
68(39)%. 



IV. eeqq CONTACT INTERACTIONS 

In this section we present the results of our global fit to the eeqq contact interactions. The 
HERA data. Tables H and mainly constrain the eeuu and eedd contact interactions, as do 
the cesium atom weak charge (|26| ) and the three asymmetries (|31|), (pS]), ( p9| ) of polarized 
electron scattering on nuclei targets. The Drell-Yan (DY) lepton-pair production data. 
Table [m] , are sensitive to eeqq and fifj-qq contact interactions mainly of the first generation 
quarks {u and d). The DY data show no indication of e-fi universality violation, hence we 
assume the e-/i universality of the contact interactions when we include this data in our 
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fit. The LEP/SLC measurements of e~^e hadrons, Table |^ and Eq. (^Hj) constrain 



eeqq contact interactions summed over five fiavors of quarks, u, d, s, c and b. Since b and 
c fiavor-tagging measurements have only limited accuracy, we need to assume quark-fiavor 
universality relations when we include the LEP/SLC data in our global fit, 

C, = C,, v:$ = Va, = V%, (55) 

where a(3 = LL, LR, RL, RR. Finally, the neutrino-nucleon scattering data ( ^9D constrain 
the neutral current v^j_v^uu and z/^z/^dd contact interactions. To include neutrino data in the 
global analysis we need to assume the electroweak gauge symmetry relations of Eq. (|^). Our 
global fits for the eeqq contact interactions are thereby organized in the increasing order of 
model dependence. 

In Table we summarize the experimental constraints obtained on the eeqq contact inter- 
actions r/^^ in the various "models" : (a) the SM result, (b) the choice r/^'^ = r/^*^: — 1-4 GeV^^ 



used in our previous qualitative demonstration of the effects of the contact interactions [|10 
(c) the fit when all 8 eeqq contact terms are allowed to vary freely, (d) a minimal choice 
that would explain the HERA and other data, (e) a vector- vector model, (f) an axial vector 
- axial vector model, and (g) the SU(12) model []TT|. The fits were performed using the 



program MINUIT [|6 



The breakdown of the minimum for each model is summarized in Table 0. The 



SM global-fit is excellent, giving x^/d.o.f. = 176.4/164. The SM gives the best of the 
CDF-DY data. Only for the HI data does the SM give a poor representation. Our previous 
choice (b) fits the HI data but does not fit ZEUS nor CDF-DY data well. For case (c) in 
which all eight parameters are allowed to vary, the minimization maintains a balance among 
the of the data sets. From Tables |V| and |V| we see that the best estimate for the LR and 
RL components are similar and the minimum for (c) and (d) also is comparable. This 
confirms that the HERA anomaly is better explained by LR and RL components, which 
have less infiuence on other experiments. In the cases (e), (f), and (g) the for the HERA 
data does not improve over (c) and (d), which indicates that the VV , AA, and the SU(12) 
chirality combinations cannot better explain the HERA anomaly without upsetting other 
experiments. Other 77^^ solutions could conceivably give a smaller for HERA but only at 
the expense of larger x^'s for the other experiments. The above fits in Table |V| and |VI| do 
not include the uN data because most of the "models" do not maintain the SU(2) invariance 
relations and (|p. 

In all the fits, we notice that rjjfpf^ and rj'^jj are in general larger than the others, even 
larger than Tjf^^ and rj^^L- This is because (i) the excess events at HERA is best explained by 
by the LR and RL chiralities (see Eq. (p^), (ii) in Drell-Yan production at the Tevatron 
uu annihilation dominates over dd, and the Drell-Yan cross sections agree with the SM 
prediction. Therefore, the minimization process increases rif^j^ and ?7|j^ to fit the HERA 
data without spoiling the SM fit to the Drell-Yan cross section. 

Table |V11| shows how the best estimates for 77^^ parameters change when data sets are 
added successively in the minimization. The most dramatic changes occur when the data 
for Drell-Yan (DY) lepton-pair production at the Tevatron are included. Note that when 



we include the z/-N data in the last column of Table |VII| the SU(2)i invariance relations of 



Eqs. (0) and (B) are assumed and then we have 7 free parameters instead of 8. Table |Vli 



shows that the HERA and DY data are not fully compatible. There is an abrupt change 
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of the best-fit 77^^ values and a sudden increase of the from the HERA data between 
the third (HERA+APV+eN) and fourth (HERA+APV+eN+DY) columns. The difference 
XsM ~ Xmin is as large as 12.8 at the third column but it reduces to 8.2 in the fourth column 
after including the CDF-DY data, which is no significant improvement over the SM since 
there are 8 eeqq contact terms allowed to vary freely. 

Figure |I] shows the fit for the HERA data and Fig. ^ shows that for the CDF-DY data 
from the models as given in Table |T|: (i) fit to HERA only (long-dashed curves), (ii) fit 
to HERA-|-APV-|-eN data (dotted curves), (iii) fit including the DY data assuming e-/i uni- 
versality (dash-dotted curves), (iv) fit including the LEP and z/N data, assuming the quark 
flavor universality and SU(2)2, invariance relations (dashed curves). The incompatibility of 
the HERA high-Q^ data and the CDF high mass DY data can be seen in these figures. Even 
allowing for general eeqq contact interactions, it is not possible to find a good fit to both 
the HERA and CDF-DY data. 

Finally, we present the result of the global fit to all low and high energy lepton-quark 
experiments by assuming the fiavor universality and the 'S\J{2)l invariance of the contact 
interactions. We then have the seven independent lepton-quark contact interactions that 
appear in the last column of Table |VI1] . The errors given for each 77^^ are correlated. The 
correlation matrix as obtained from MINUIT PBI is as follows: 



/ 


1.00 


-0.65 


0.78 


-0.52 


0.91 


0.14 


-0.13 \ 




-0.65 


1.00 


-0.76 


0.69 


-0.38 


-0.31 


0.47 




0.78 


-0.76 


1.00 


-0.71 


0.58 


0.36 


-0.49 




-0.52 


0.69 


-0.71 


1.00 


-0.43 


0.22 


0.75 




0.91 


-0.38 


0.58 


-0.43 


1.00 


-0.18 


-0.10 




0.14 


-0.31 


0.36 


0.22 


-0.18 


1.00 


0.35 


V 


-0.13 


0.47 


-0.49 


0.75 


-0.10 


0.35 


1.00 / 



(56) 



where the row and column indices are labeled in the order: t]!*^, t^I'^, 77^')^, r;^, ?7|5^, T^^'fj, 
Tj^YiR- The covariance matrix C can be formed by Cij = aiajpij, where ai is the error of the 
i-th parameter. The eigenvectors of the covariance matrix are constrained as: 



.7377/21 + 
.2167/21 - 
.0357/21 - 
.4297/21 + 
.0637/21 - 
•268r/21 - 



•3197/21 - 

•13lCi? + 
.5457/21 - 
•6 787/21 - 
•0417721 + 
.0987/21 + 
•336r/21 + 



.1267/- + 
.2327/- + 
,394r/|jl - 
.0967/- + 
.1877/- + 
.5907/- - 



•0807/- - 
.6827/- + 
•4737/le"^ - 
•0707/- + 
•1587/^"« - 
•434//- + 



•26l7/2i + 
.6307/21 - 
.1537/21 - 
.1997/21 - 
.6387/21 - 
.0537/21 + 
•249r/2l + 



•3217/21 - 
.1337/21 + 

•014r?2l 
•0477/21 + 
•2307/21 + 
.9047/21 + 
•0 797/21 - 



•329<^ 
•0557/^1 
•3387/^1 
•3437/1,1, 
•5837/^1 
•3267/^1, 
•459<« 



0.034 + 0.031, (57a) 

0.014 + 0.063, (57b) 

-0.92 + 0.29, (57c) 

-0.45 + 0.37, (57d) 

0.33 + 0.82, (57e) 

1.04+ 1.46, (57f) 

-0.36 + 1.50. (57g) 



The strongest constraint ( p7a| ) is essentially due to the atomic parity violation experiment 
and hence it agrees roughly with the result of the corresponding fit to the low energy 
electroweak data only The rest of the constraints show significant improvements over 
those obtained from low energy experiments, even though the analysis of Ref. [38| made 
a more restrictive assumption that r/21 = //21- The combined effect of low-energy and 
high-energy data is that irrespective of the fiavor and chirality combination no eeqq contact 
interactions can be significantly larger than about 1 TeV~^. 
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Throughout, we have parameterized the eeqq contact interactions via the rf^f^ which, to a 
reasonable approximation, are normally distributed. Constraints on contact terms are often 
quoted in terms of 95% CL bounds on the contact scale A^'^^ which is related to the r^'s via 

where e = ±l,a,/5 = L, i?. In order to allow an easy comparison we also derive 95% CL 
bounds on the A^"^. The procedures for obtaining these limits is as follows: 

1. Using MINUIT we obtain the best estimate of a parameter r^^^ (others set to zero) by 
minimizing the fo^' the available data. When obtaining t]^^ = ± a, we assume 
that the probability density function P{ri) follows a normal distribution with mean ii 
and variance cr^: 

P(r,) = -±= e-«^)" . (59) 
av ATI 



Possible non-Gaussian behaviour is taken into account by replacing (2-Ji)2 -\^^ x^iv) " 
Xmin argument of the exponential. 

We determine "+" and "—"-sided limits of A^'^^. The "+"-sided limit means that the 
physically allowed region of rf^^ is rf^^ > 0. To obtain a 95% CL limit on rj'^^ we need 
to find the value, 77+, above which 5% of the total probability for positive t] is located, 

r° PMdn 
/o P{v)dv 



This 1]+ is related to A+ by 77+ = 47r/(A 



The "—"-sided limit means that the physically allowed region of rj^p is rj'^n < 0. To 



obtain a 95% CL limit on r/^^ we need to find the value, below which 5% of the 
total probability for negative rj is located, 

0.05 = )'i ' . (61) 

This ?7_ is related to A_ by r]^ = —An/ (A_)^. 

The limits are tabulated in Table |V111| for individual couplings of a fixed lepton and quark 
chirality. The SU(2) relations of Eqs. (^ and (^ are imposed and only one contact term at 
a time is assumed to be different from zero. Combinations of contact interactions which are 



parity conserving or blind to atomic parity violation constraints are considered in Tables 
and ^ 
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V. CONCLUSIONS 



In summary, we have made a general investigation to determine whether the excess of 
events observed at high-Q^ in the HERA experiments can be understood in terms of contact 
interactions. All available low and high energy data relevant to eeqq contact terms have 
been included in a global fit. Our findings can be summarized as follows: 

(i) If the HERA high events are due to eeqq contact interactions, rjf^p^ and rj'j^]^ terms are 
the most effective in explaining the data. 

(ii) The CDF data on high mass Drell-Yan lepton-pair production, however, strongly restrict 
the eeuu contact interactions irrespective of their chirality structure. 

(iii) Once the high-energy data and the low-energy electroweak data are combined, all eeqq 
contact interactions are strongly constrained and the possibility of explaining the HERA 
high events with contact interactions is limited. 

(iv) The SM gives a reasonably good overall description of all the electron-hadron data, 
including the HERA high events. 

(v) The limits obtained on the contact scales are more restrictive than those obtained pre- 
viously p5| . 

ACKNOWLEDGMENTS 

K.C would like to thank S. Godfrey and K. McFarland for discussions. We thank 
Dieter Haidt for discussions of error correlations in the HERA data. This research was 
supported in part by the JSPS-NSF Joint Research Project, in part by the U.S. Department 
of Energy under Grant Nos. DE-FG03-93ER40757 and DE-FG02-95ER40896 and in part 
by the University of Wisconsin Research Committee with funds granted by the Wisconsin 
Alumni Research Foundation. 



17 



REFERENCES 



[1] 
[2] 
[3] 

[4] 
[5] 



The HI Collaboration, C. Adloff et al., Z. Phys. C74, 191 (1997). 

The ZEUS Collaboration, J. Breitweg et al, Z. Phys. C74, 207 (1997). 

G. Aharelli, J. Ellis, G.F. Guidice, S. Lola, and M.L. Mangano, report CERN-TH/97- 

040 ||hep-ph/970327'qi . 

K.S. Babu, C. Kolda, J. March-Russell, and F. Wilczek, Phys. Lett. B402, 367 (1997). 
J. Bliimlein, Z. Phys. C74, 605 (1997); J. Kalinowski, R. Riickl, H. Spiesberger, and 
P.M. Zerwas, Z. Ph ys. C74, 595 (19 97); M. Suzuki, |hep-ph/97033T6| ; G.K. Leontaris 
and J.D. Vergados, |hep-ph/9703338| ; Z. Kunszt and W.J. Stirling, Z. Phys. C75, 453 
(1997); T. Plehn et al, Z. Phys. C74, 611 (1997); D. Friberg, E. Norrbin, and T 
Sjostrand, Phys. Lett. B403, 329 (1997); I. Montvay, Phys. Lett. B407, 22 (1997); S.F 
King and G.K. Leontaris, [hep-ph/9704336| ; J. Elwood and A. Faraggi, |hep-ph/ 9704363 
B. Dutta, R. Mohapatra, and S. Nandi, |hep-ph/9704428| 



S. Jadach, W. Placzek, and B 



Ward, |hep-ph/ 97053^ ; K.S. Babu, C. Kolda, and J. March-Russell, |hep-ph/9705411 



J. Ellis, S. Lola, and K. Sridhar, |hep-ph/97054T6| ; S. Jadach et al., |hep-ph/9705429| ; A 
Blumhofer and B. Lampe, |hep-ph/ 970615^ ; E. Keith and E. Ma, |hep-ph/9707'2T^ ; N.G 



[6] 
[7] 



Deshpande and B. Dutta, |hep-ph/9707274 



[9] 
[10] 

[11] 
[12] 

[13] 

[14] 

[15] 



[16] 

[17] 
[18] 

[19] 

[20] 



T. Kon and T. Kobayashi, Phys. Lett. B270, 91 (1991). 

D. Choudhury and S. Raychaudhuri, Phys. Lett. B401, 54 (1997); H. Dreiner and P. 
Morawitz, |hep-ph/9703279| ; J. Kalinowski, R. Riickl, H. Spiesberger, and P.M. Zerwas, 
Phys. Lett. B406, 314 (1997); T. Kon and T. Kobayashi, |hep-ph/ 9704221] ; R. Bar- 
bieri, A. Strumia, and Z. Berezhiani, [lep-ph / 9704275] ; G.F. Giudice and R. Rattazzi, 
Phys. Lett. B406, 321 (1997); A.S. Belyaev and A.V. Gladyshev, |hep-ph/9704515| ; J. 
Kalinowski, |hep-ph/9706203| ; J.E. Kim and P. Ko, |hep-ph/9706387| ; A.S. Joshipura, 
V. Ravindran, and S.K. Vempati, |hep-ph / 9706482] ; M. Frank and H. Hamidian, ]hep-] 
ph/97065iq ; M. Guchait and P .P. Roy, jhep-ph/ 9707275] ; T. Kon, T. Matsushita, and 
T. Kobayashi, ]hep-ph/9707355 . 

S. Adler, ]hep-ph / 9702378] ; K. Akama, K. Katsuura, and H. Terazawa, lhep-ph/9704327 . 

C. J. C. Purges and H. J. Schnitzer, Phys. Lett. 134B, 329 (1984). 

V. Barger, K. Cheung, K. Hagiwara, and D. Zeppenfeld, Phys. Lett. B404, 147 (1997). 

A. Nelson, Phys. Rev. Lett. 78, 4159 (1997). 

N. Bartolomeo and M. Fabbrichesi, Phys. Lett. B406, 237 (1997). 
W. Buchmiiller and D. Wyler, ]hep-ph/97043T7 . 

N.G. Deshpande, B. Dutta, and Xiao-Gang He, preprint OITS-630 |]hep-ph/9705236| |. 
M.C. Gonzalez- Garcia and S.F. Novaes, ]hep-ph / 9703346] ; S. Godfrey, Mod. Phys. Lett. 
A12, 1859 (1997); L. Giusti and A. Strumia, ^iep-ph/9706298] ; Z. Cao, X.-G. He, and 

B. McKellar, ]hep-ph/9707227] ; F. Cornet and J. Rico, ]hep-ph/ 9707229 



S. Kuhlmann, H.L. Lai, and W.K. Tung, ]hep-ph/97043381 ; K.S. Babu, C. Kolda, and J. 
March-Russell, ]hep-ph/9705"399] ; J.F. Gunion and R. Vogt, ]hep-ph/9706^ . 
M. Drees, Phys. Lett. B403, 353 (1997); U. Bassler and G. Bernardi, ]hep-ex/9707024 . 
H. Kambara for the CDF Collaboration, talk given at Hadron Collider Physics XII, 
Stony Brook, June 5-11, 1997. 

D. Norman for the DO Collaboration, talk given at Hadron Collider Physics XII, Stony 

Brook, June 5-11, 1997. 

C.S. Wood et al.. Science 275, 1759 (1997). 



18 



[21] C.Y. Prescott et al., Phys. Lett. B84, 524 (1979). 

[22] W. Heil et al., Nucl. Phys. B327, 1 (1989). 

[23] P.A. Souder et al., Phys. Rev. Lett. 65, 694 (1990). 

[24] A. Bodek for the CDF Collaboration, preprint Fermilab-Conf-96/341-E (1996); CDF 
Collaboration, F. Abe et al., Fermilab-PUB-97/171-E. In our analysis, we used the data 
from Fig. 1 of the latter paper. 

[25] LEP Collaborations and SLD Collaboration, "A Combination of Preliminary Elec- 
troweak Measurements and Constrains on the Standard Model" , prepared from contri- 
butions to the 28th International Conference on High Energy Physics, Warsaw, Poland, 
CERN-PPE/96-183 (Dec. 1996). 

[26] OPAL Collaboration, G. Alexander et al., Phys. Lett. B376, 232 (1996); ibid B387, 
432 (1996); ibid B391, 221 (1996). 

[27] L3 Collaboration, Phys. Lett. B370, 195 (1996); CERN-PPE/97-52, L3 preprint 117 
(May 1997). 

[28] ALEPH Collaboration, Phys. Lett. B378, 373 (1996). 

[29] P. Langacker and J. Erler, presented at the Ringberg Workshop on the Higgs Puzzle, 

Ringberg, Germany, 12/96, |Ei^h/9703428| . 
[30] K.S. McFarland et al. (CCFR), FNAL-Pub-97/OOl-E, |hep-ex/970lOTq . 
[31] E. Eichten, K. Lane, and M. Peskin, Phys. Rev. Lett. 50, 811 (1982). 
[32] R.J. Cashmore et al., Phys. Rev. 122, 275 (1985); R. Riickl, Phys. Lett. B129, 363 

(1983); Nucl. Phys. B234, 91 (1984). 
[33] P. Chiappetta and J.-M. Virey, Phys. Lett. B389, 89 (1996). 

[34] C. J. C. Burges and H. J. Schnitzer, Nucl. Phys. B228, 464 (1983); W. Buchmliller and 

D. Wyler, Nucl. Phys. B268, 621 (1986). 
[35] Review of Particle Properties, Phys. Rev. D54, 1 (1996). 

[36] W. Buchmiiller, R. Riickl, and D. Wyler, Phys. Lett. B191, 442 (1987); J. Hewett and 

T. Rizzo, [hep-ph/9703337| , to be published in Phys. Rev. D. 
[37] P. Langacker, M. Luo, and A. Mann, Rev. Mod. Phys. 64, 87 (1992). 
[38] G.-C. Cho, K. Hagiwara and S. Matsumoto, |hep-ph/970733^ . 

[39] K. Hagiwara, D. Haidt, C.S. Kim and S. Matsumoto, Z. Phys. C64, 559 (1994); Errata 

C 68, 352 (1995). 
[40] D.H. Beck, Phys. Rev. D39, 3248 (1987). 
[41] H. Burkhardt and B. Pietrzyk, Phys. Lett. B356, 398 (1995). 

[42] V. Barger and R.J.N. Phillips, Collider Physics, (Addison- Wesley, Reading MA, up- 
dated edition 1996). 

[43] S.G. Gorishny, A. Kataev, and S.A. Larin, Phys. Lett. B259, 114 (1991); L.R. Surgu- 

ladze and M.A. Samuel, Phys. Rev. Lett. 66, 560 (1991). 
[44] G. Fogli and D. Haidt, Z. Phys. C40, 379 (1988). 

[45] A. Blondel et al. (CDHS), Z. Phys. C45, 361 (1990); J.V. Allaby et al. (CHARM), Z. 

Phys. C36, 611 (1987); D. Allasia et al. (BEBC), Nucl. Phys. B307, 1 (1988); C.G. 

Arroyo et al. (CCFR), Phys. Rev. Lett. 72, 3452 (1994). 
[46] CERN Program Library entry D506. 



19 




FIG. 1. Integrated cross sections for e'^p — > e^X versus a minimum for the SM (solid 
curve) and for four choices of contact interactions: best-fit for the HI and ZEUS data (long-dashed 
curved), best-fit including the low energy electron-quark coupling measurements (dotted curve), 
best-fit further including the Tevatron lepton-pair production data (dot-dashed curve), and the 
best-fit including also the LEP data and the neutrino-scattering data (dashed curve); see Table 



VI] . The data points correspond to the number of events observed by HI |1[ and ZEUS M, and 



have been converted to cross sections by us, assuming constant detection efficiencies of 80% and 
81.5%, respectively. Our fits are performed directly on the published data im2| as summarized in 
Tables | and 0. 



20 




FIG. 2. The Drell-Yan cross section pp — > or e~^e~)X at the Tevatron {^/s = 1.8 TeV) 

for the SM (soUd curve) and the four choices of contact interactions as in Fig. |l[ Prehminary CDF 
data |25| are shown separately for e+e" (soHd triangle) and (open square) pair; see Table III. 

The cross section is averaged over y as ^ /j*;^ dy jj[/^y ■ 
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TABLE I. The expected and observed numbers of events in each bin of the x — y plane from 
ZEUS. The expected number is hsted above the observed number. 
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TABLE II. The expected and observed numbers of events in each bin of the x — y plane from 
HI. The expected number is listed above the observed number. 
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0.3592 

0.3791 


0.3794 

0.1002 


0.4002 

0.1216 


0.4216 

0.1135 


0.4435 

0.95 


0.5 < y < 0.9 


73.59 
76 


2.2 

3 


0.27 



0.21 



0.15 
1 


0.15 



0.14 
1 


0.39 
3 


0.4 < 7/ < 0.5 


56.28 
59 


2.28 
2 


0.25 



0.21 



0.19 



0.07 
1 


0.08 
1 


0.24 



0.3 < y < 0.4 


65.13 

80 


3.77 

3 


0.40 
1 


0.42 
1 


0.36 



0.24 



0.19 




0.69 




0.2 < y < 0.3 


70.9 
64 


9.93 
12 


1.05 
2 


0.91 



0.75 
1 


0.61 
1 


1.45 



1.4 
1 



22 



TABLE III. Preliminary CDF Drell-Yan Data extracted from their figures. Both e and 
[i^ iT' samples are shown. The data is the y-averaged integrated cross section ^ /^^ f^y-^'^ " 



' dMdy- 



Mee(GeV) 


i/-i^l/dS;(pb/GeV) 


M^^(GeV) 


|/-i^^JS;(pb/GeV) 


54.63 


0.142 ±Hi 


54.63 


0.118 ±0.019 


64.90 


+0.017 

0.113 -0.016 


64.90 


n nno +0-015 
0.098 -0.014 


73.82 


4-0 01 fi 

0.130 -om7 


73.82 


0.113 -oill 


81.86 


1 Q Q29 

0.358 -0^035 


81.86 


0.329 -0.030 


87.66 


_i_n 1 fift 

2.537 -oiti 


87.66 


t n 1 KQ 

2.162 ioii 


91.68 


^ „ „„„ +0.861 
10.779 -0.668 


91.68 


p. +0.906 
11.345 -0.839 


97.92 


0.462 -0.045 


97.92 


+0.044 

0.502 -0.040 


105.96 


+0.022 

0.122 -0.020 


107.74 


, +0.020 

0.148 -0.021 


114.88 


n n/lov +0-0104 
() ()4z / —0 nio^ 


116.67 


n nroK +0 0125 
( ( h,-)h — n 01 31 


134.97 


0.0209 ±0.0041 


138.98 


0.0214 IHSi 


174.69 


0.00352 ±0.00127 


178.71 


0.00426 


225.11 


0.00234 ±0.00105 


228.69 


n nnintx +0-00099 
0.00 195 -0.00097 


275.11 


(9.41 ±6.70) X 10-^ 


279.13 


(9.91 ± 7.05) X 10-^ 


350.09 


(^2.22 -lil^ X 10-^ 


354.11 


(2.22 111) X 10-4 


450.06 


(2.27 ±2.27) X 10-"^ 


454.08 


(0.0 t'o^o) X 10-4 


550.04 


(^0.0 t'o'o') X 10-^ 


554.06 


(0.0 t'oL^j X 10-4 



TABLE IV. Table showing the LEPl and SLD data, and LEP1.5 and LEP2 data that are 



relevant to our analysis |25-2£] 



EcM 


0"had 


SM 
'^had 


LEPl, SLD ■.^/^ = Mz 


41.508 ± 0.056 nb 


41.465 nb 


OPAL v 


= 130.26 GeV 


66 ± 5 ± 3 pb 


78 pb 


OPAL v 


= 136.23 GeV 


60 ± 5 ± 2 pb 


63 pb 


OPAL v 


= 140 GeV 


50 ± 36 ± 2 pb 


56 pb 


L3 


V~s = 


130.3 GeV 


81.8 ±6.4 pb 


78 pb 


L3 


V~s = 


136.3 GeV 


70.5 ± 6.2 pb 


64 pb 


L3 


V~s = 


140.2 GeV 


67 ± 47 pb 


56 pb 


ALEPH : 


y/s = 130 GeV 


74.2 ± 5.2 ± 3.3 pb 


76.9 pb 


ALEPH : 


v/i = 136 GeV 


57.4 ± 4.5 ± 1.8 pb 


62.5 pb 


OPAL v 


^ = 161.3 GeV 


35.3 ± 2.0 ± 0.7 pb 


33.2 pb 


L3 


V~s = 


161.3 GeV 


37.3 ±2.2 pb 


34.9 pb 


L3 


V~s = 


170.3 GeV 


39.5 ± 7.5 pb 


29.8 pb 


L3 


V~s = 


172.3 GeV 


28.2 ± 2.2 pb 


28.9 pb 
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TABLE V. The best estimate of the contact interaction parameters ry^^ and the corresponding 
Xmin i*^) eight ?7^^ being free, and (d)-(f) various chirahty combinations. The t^ie 
Standard Model and the "eye-ball" solution from our previous paper are shown in (a) and (b), 
respectively. The data sets included are: HI, ZEUS, APV, e-N, DY, and LEP. The number of 
degree of freedom for the Standard Model is 164. 



Chirality Combination 



Best estimate (unit=TeV 



A mm 



[a) Standard Model 



176.4/164 d.o.f. 



(b) r?r^ = 
others=0 



1.4TeV" 



0.16 -0.48,7/1"^ 

: 0.80 iai, 7]I^L 
-1.01, r/lf^ 



187.1/163 d.o.f. 



'c) All free 



Vll 

vfR 



-0.001 foil 
0.19 trrf. 



0.31 



0.62 nil, 
1.66 tin vltL = 1-95 tin 



167.2/156 d.o.f. 



(d) r?!"^ = r]l% r/lf^ = rjl% 
others=0 



Vrl 

ed 



vTr 

ed 

Wr 



0.51 
1.97 



hO.26 
-0.27, 

h0.50 
-0.70, 



170.0/ 162 d.o.f. 



(e) Ti^y,r]^v 



Vvv 



1-0.20 



0.47 -6:23, 7]^^y = 1.10 



-0.28 
-0.38 



172.5/162 d.o.f. 



Vaa 



0.17 



-0.36 -5:i4, r]f^ 



-0.36 -0.46 



173.0 /162 d.o.f. 



:g) su(i2) 



Vll 

eu 

Vrl 

ed 

Vll 
V'A 



-Vlr ■ 

-vWr 

ed 

-Vlr ■ 
-vtR 



n on +0-68 
-0.39 -0.31, 



0.29 



[-0.48 
-0.68, 



-0.81 ioJe, 
-0.11 



172.7/160 d.o.f. 



TABLE VL Breakdown of the minimum by experiment for the models in Table 





HI 


ZEUS 


APV 


eN 


CDF-DY 


LEP 


# Data points 


32 


81 


1 


4 


32 


14 


(a) SM 


33.9 


59.9 


1.2 


1.8 


62.2 


17.4 


(b) vTr = VtL = 1.4TeV-2, 
others=0 


29.0 


66.5 


1.2 


1.8 


70.9 


17.7 


(c) All free 


29.9 


58.1 


0.0 


0.5 


62.6 


16.1 


(d) r?|e'i = r?r^ = 0.51, 
rif^ = r]fj^ = 1.97 


30.0 


58.2 


1.2 


1.8 


62.7 


16.1 


(e) r/PV = 0.47, 7/1.^ = 1-10 


30.4 


58.7 


1.2 


1.8 


63.1 


17.4 


(f) r/2^^ = -0.36,<^ = -0.36 


31.9 


58.8 


1.2 


1.8 


62.2 


17.1 


(g) vTl = -vTr = -0.39, 
^ITl = = 0.29, 
Vli = -vtR = -0.81, 

= -vA = -0.11 


32.0 


58.8 


1.2 


1.5 


62.3 


17.0 
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TABLE VII. The best estimate of the rf^^ parameters when various data sets are added 
successively. In the last column when the f-N data are included the t/^^ are given in terms of rf^^ 
by Eq. @, we assume ryljj^ = r]'^^^ in the last column. 





HERA only 


HERA-fAPV 


HERA-fAPV 


HERA-FAPV 


HERA-FAPV 


HERA-hDY-FAPV 








+eN 


-heN-FDY 


-FeN+DY-FLEP 


-FeN-FLEP-hi/N 






omtlil 


0.49til^ 


0.16-0.53 


0.16-0.48 


-0.082^0:40 


ILK 


-4.93io:96 


-4.90i|72 


-2.24lt:^^ 


o.88lg:i 




o.86lg:i 




-U. 14—5.07 


-U.oy— 4.80 


-4.ZU— 1.08 


U.zi— 0.78 


u.iy— 0.71 


U.oo— 0.91 


^RR 


+4.42 
1.28-6.31 


+4.45 
1.03-5.90 


. +1.73 
1.46-3.73 


+0.65 
0.011-0.61 


+0.63 
-0.0010-0.55 


^.+0.74 

-0.074-0.57 




-2.51 -8.24 


-i.88i?:59 


-4.o3iI:oi 


o.7iil:i 


o.3il?:o? 


0.015^0:1 


Wr 


-1.42-4.87 


-1.19-4.95 


-1.97-3.77 


1.15-1.99 


1.66-1.66 


0.88-1.82 


Vrl 


-2.8oi5:?6 


-3.34ll:i 




i.sollJ? 


1.95lli^ 


= Vrl 


'Hri 


-3.16 -9.18 


-4.13l|i 




i.oiil:P 




o.84l?:?f 


HERA 


83.4 


83.4 


83.6 


88.1 


88.0 


88.6 


APV 




0.0 


0.0 


0.001 


0.0 


0.001 


eN 






0.38 


0.47 


0.47 


0.67 


DY 








62.3 


62.6 


62.2 


LEP 










16.1 


17.1 


i/N 












0.007 


Total 


83.4 


83.4 


84.0 


150.8 


167.2 


168.7 


SM 


93.8 


95.0 


96.8 


159.0 


176.4 


176.4 


SM d.o.f. 


113 


114 


118 


150 


164 


165 



TABLE VIII. The best estimate on -rf^^ and the 95% CL limits on the contact interaction 
scale A^'^, where rf^^ = 47re/(A^'^^)^. When one of the r/'s is considered the others are set to zero. 
SU(2) relations are assumed and z^N data are included. 



Chirality [q) 


r/ (TeV-2) 


95% CL Limits 
A+ (TeV) A_ (TeV) 


LL(n) 


0.046 ± 0.057 


9.3 


12.0 


LR(ii) 


0.11 ±0.079 


7.2 


11.3 


RL(u) 


-0.041 ± 0.038 


15.4 


10.9 


RR(u) 


-0.11 ±0.080 


11.1 


7.3 


LL(d) 


0.055 ± 0.061 


8.8 


11.9 


LR(d) 


0.076 ± 0.072 


7.9 


11.2 


RR((i) 


-0.074 ± 0.072 


11.2 


8.0 
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TABLE IX. The best estimate on r/'^'^ for the minimal setting, VV,AA, and SU(12), and the 
corresponding 95% CL Hmits on the contact interaction scale A, where rj = 47re/(Ag)^. When one 
of the rj's is considered the others are set to zero. Here we do not use SU(2) relations nor do we 
include the uN data. 







yo/o UL; Liimits 




Chirality (q) 


ri (TeV-2) 


A+ (TeV) 


A_ (TeV) 


Vlr = Vrl 


0.58 t'oE 


3.6 


5.8 


Vlr = Vrl 


2.07 t'oil 


2.1 


2.6 


V^v 


0.024loii 


6.2 


7.5 


vfv 


0.29tg:P 


3.8 


2.6 


Vaa 


-0.28lgi! 


6.5 


5.1 


Vaa 


0.21 t'oM 


4.3 


5.1 


Vll = -Vlr 


-QA2t'ol 


4.2 


4.1 


'IrL — 'IRR 


0.51-0.26 


3.9 


3.9 


vii = -vtR 


o.3itg:lg 


3.4 


2.2 


Vrl = -Vrr 


—0.51-0.50 


2.3 


3.1 



TABLE X. Same as Table II but with a further condition: 77^" = 77^". Here q = u = d. 







9 


5% CL Limits | 


Chirality (g) 


rj (TeV-2) 


A+ (TeV) 


A_ (TeV) 


eq eq 

Vlr = Vrl 


0.69 


3.5 


5.3 


eq 

Vvv 


o.iitg:^ 


4.1 


6.9 


eq 

Vaa 




4.4 


4.7 


eq eq 

Vll = -Vlr 

eq eq 
Vrl — Vrr 


-0.54-0.20 

0.58-0.26 


3.5 
3.8 


3.8 
3.4 
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